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ABSTRACT 


Simultaneous measurements of sea surface elevation and onshore and 
alongshore water particle velocities were measured at three locations 
within the surf zone using two capacitance type pentrating wave staffs 
and three two-component electromagnetic flow meters. The probability 
density functions, pdf, for the sea surface elevation were always highly 
positively skewed, whereas the pdf's for the velocities were both nega- 
tively and positively skewed. Mean values of the onshore and alongshore 
components of flow reflected the influence of a rip current frequently 
observed just south of the instrument locations. Strong harmonics in 
the spectra of sea surface fluctuations and particle velocities infer 
nonlinear conditions. Coherence values between waves and onshore flow 
were high, ranging above 0.9. The coherence between waves and onshore 
flow was used to separate the turbulence and wave-induced velocity 
components. Over the range of collapsing to spilling breakers a rea- 
Sonable value for the ratio of turbulent to wave-induced velocity was 
determined to be approximately 0.75. Saturation regions were found in 
the wave and velocity energy-density spectra at higher frequencies as 


evidenced by -5 and -3 slopes, respectively. 
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I. INTRODUCTION 


С TOPICAL PERSPECTIVE 

The surf zone is an area bounded on the seaward side by the point 
where waves first begin to break, and on the landward side by the point 
of maximum run-up on the beach slope. Although beach profiles and waves 
within the surf zone have been studied by many, the kinematics of wave 
forms and water particle velocities in the zone have remained something 
of a mystery. 

The problems encountered have been practical as well as theoretical. 
Wave theories developed for deep water waves do not correctly characterize 
the motion which occurs as the waves break. Direct measurements of break- 
ing waves appear to be the most viable means of approaching the problem. 
However, this has been hampered by the surf zone's hostile environment, 
which is difficult to reproduce in the laboratory, and by inadequate 
instrument design. Improvements in the latter area has resulted in sturdy, 
sensitive, but expensive measuring devices which have rapid response times. 

The earliest study of significant importance was conducted by 
Iversen (1953) in which he used photographic techniques to obtain a 
Lagrangian description of water particle motion under breaking waves. 
The laboratory channel limited the wave type to plunging and surging 
breakers. 

Inman (1956) was one of the first to make sophisticated field 
measurements, measuring the drag force on a cantilevered sphere to infer 


water particle motion. Miller and Zeigler (1964) measured in situ 








particle motions using acoustic and electromagnetic current meters, 

then compared their findings with higher order wave theory and found 

some qualitative agreement. Walker (1969) made similar measurements 
using propeller type flow meters. Meadows (1976) used ducted impeller 
flow meters at equally spaced vertical positions in order to measure 
longshore components of flow. Huntley (1976) utilized a single two- 
component electromagnetic flow meter in an effort to obtain a value 

of the friction coefficient. Bowen and Huntley (1974) made measure- 
ments of the nearshore velocity fields using as many as three two- 
component electromagnetic flow meters. Führböter and Büsching (1974) 
measured simultaneous orbital velocities and water levels using a two- 
component current meter and two pressure type wave meters. Thornton 
(1969), Steer (1972), Thornton and Richardson (1973), Bub (1974) and 
Galvin (1975) used pressure meters, capacitace wave gauges and electro- 
magnetic current meters tocompute particle velocities and surface profiles 
within the surf zone. The work presented herein is an extension of these 


latter studies. 


В. OBJECTIVE 

The objective of this research is to study the kinematics of water 
particles within breaking waves and within the surf zone. In particular, 
the experiment was performed on a flat beach in order to measure the 
characteristics of spilling breakers. Simultaneous measurements were 
made of tne instaneous sea surface elevation and two orthogonal water 
particle velocities, Onshore and alongshore, at three fixed locations 


in the surf zone. Estimates of the probability density functions and 
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power spectra of the wave heights and particle velocities are made. The 
suitability of using linear wave theory as a Spectral transfer function 
in determining velocity spectral components from the power spectrum of 
the waves is measured. Computed velocity spectra are compared with 
actually measured velocity spectra. Coherence between waves and onshore 
flow is used to separate the turbulence and wave-induced velocity com- 
ponents in order to determine the ratio of turbulent to wave-induced 
velocity. Spectral estimates are analyzed to ascertain the slope of 


the saturation region at higher frequencies. 
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11. MEASUREMENTS 


A. EXPERIMENT 

The experimental site is just north of La Jolla, California. Data 
presented in this study was taken as part of a larger scale experiment 
conducted from 15 February to30 March 1977 in the Southern California 
bight area. The over-all purpose of the experiment was to test and 
evaluate the instrumentation package which was selected for use aboard 
SEASAT. The apparatus, which will measure waves, winds, and sea surface 
temperature, was being flown on board NASA airplanes. A major ground 
truth program was concurrently being conducted at Torrey Pines Beach. 

Measurements of waves, currents, and sediment transport at Torrey 
Pines Beacn were made using the following equipment. A five element 
linear array of pressure sensors was positioned offshore along the ten 
meter contour to measure wave energy and direction. Located in the surf 
zone inshore of the array were three support towers in a line perpendi- 
cular to the beach. A surface piercing capacitance wave gauge and a two- 
component electromagnetic flow meter were mounted on each tower. A 
nephelometer was installed on the middle support tower. Te the north 
of the three towers was an array of four flow meters which was used to 
delineate the areal and temporal variation in alongshore currents. 
Figure 1 shows instrument locations. Resistance-wire run-up meters were 
installed in the swash zone at four locations. A11 instrumentation 
was battery operated. 

The offshore pressure sensor outputs were telemetered directly to 


the shore processes lab at Scripps Institute of Oceanography, which is 
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located one mile to the south. All data from the other instruments 

were cabled to one of two transmitting terminals on the beach where it was 
then telemetered to the shore processes lab and recorded. All the equip- 
ment except the nephelometer were in position prior to high tide on 9 
March. The nepholometer was installed on 14 and 15 March. Data was 
saumu lated during high tide on 9, 10, 11, 16, 17, 18, 19, 21 and 23 
March. 

Several other smaller scale experiments were conducted at the same 
time. At an elevation of approximately 300 feet on a cliff overlooking 
Torrey Pines Beach a radar installation was established. The radar 
images of waves were then photographed and analyzed. 

Sediment transport experiments were conducted on 11, 21, апа 23 
March. Bed load transport was measured using the methods described by 
Inman and Komar (1970) to trace the movement of fluorescent dyed sand. 


Suspended sediments were measured in situ by swimmers uSing a mechanical 





water sampling devide. These samples were taken along with nepholometer 
readings in order to compare methods. During these sediment transport 
investigations, Lagrangian floats made from wine bottles weighted by 


Sand were used to determine average long-shore current speeds. 


B. INSTRUMENTATION. 

The manner in which waves break depends very much on the character- 
istics of the beach and near shore bottom slope. (Table I) Spilling breakers 
Occurred most frequently at the Torrey Pines Beach site. Figure 2 shows 
a typical beach profile and the location of the instruments which provided 


the data for this paper. 
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TABLE I. BEACH AND WAVE CHARACTERISTICS 


LOCATION 

DATE 

BEACH SLOPE 

SAND TYPE 

SAND MEDIAN GRAIN SIZE 
BREAKER TYPE 

WAVE PERIOD 

WAVE HEIGHT 


TORREY PINES BEACH 

9-23 MARCH 1977 

FIGURE 2 and APPENDIX A 
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|. Wave Gauge 


The wave gauges were of the capacitance type fashioned from 
3/8 inch outside diameter stainless steel rod. The rod was tightly 
covered with 1/16 inch wall thickness polypropylene tubing. This li- 
near, highly sensitive instrument has proven to be sufficiently sturdy 
to withstand the rigors of the unfriendly surf zone. The gauges operate 
on the principle that a change in the plate dimension of the capacitor 
changes its capacitance and consequently the circuitry output voltage. 
In these gauges, the insulated steel rod and sea water act as the plates 
and the insulation functions as the dielectric. As the water level fluc- 
tuated, the capacitance changed. Fluctuations were sensed by a transis- 
torized circuit powered from the beach. The circuit was designed by 
McGoldrick (1969). The electronics packaged was housed in a watertight 
brass case which was mounted on the tower. This allowed the connecting 
leads to be relatively short, hence minimizing wire-to-wire capacitance. 
All gauges were statically calibrated in the laboratory prior to the ex- 
periment. Accuracy was estimated to Бе *.005т. The calibration plots 


are shown in Appendix B. 
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2. Flow Meter 
The flow meters were Marsh-McBirney model 511 electromagnetic 
water current meters. Operation of the meter 15 based on Faraday's prin- 
ciple of electromagnetic induction. The meter measures water particle 
velocity in two orthogonal directions through a range of zero to three 
m/sec with a maximum output error of two percent of full scale reading. 
The sensor has a variable time constant, a setting of 0.2 sec. was used. 

The flow meters were dynamically calibrated using the method of 
Thornton апа Кгароһ1 (1974). This procedure involves utilizing an os- 
cillating platform driven by a variable speed motor geared to an eccen- 
tric throw arm. This technique was employed in order to determine the 
meter's characteristics. Measurement accuracy was determined to be 
t.005 m/sec. The only difficulty encountered in their utilization oc- 
curred when Flow Meter £1 was buried on 19 March due to the heavy sedi- 
ment transport. 

The measuring instruments were attached to a 3.6m high tower con- 
structed of steel pipe with an outside diameter of 6.3cm.A 0.5m diameter 
base plate was placed about 0.6m from the bottom of the tower. This 
configuration allowed the tower to be sunk approximately 0.7m into the 
sand. The towers were also supported by steel guy wires fastened to 
blade anchors driven into the sand. Thus, tower movement and vibration 
were negligible. 

The towers were placed on a line perpendicular to the shore at 
specified distances, and were erected during low tide when the beach 


was more accessible. As indicated earlier, measurements were taken at 


18 





high tide. Instruments were arranged so as to be in the same vertical 
plane. A carpenter's level was used to establish axis alignment with an 


estimated error of 22 degrees. A typical tower and instrument arrangement 


is shown in Figure 3. 
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111. ANALYSIS OF DATA 


The data were digitized at a rate of four samples per second, 
corresponding to a sample interval of 0.25 seconds. This resulted in a 
Nyquist frequence of 2.0 Hz. The sampling rate was considered suffi- 
ciently high enough to avoid aliasing of energy into the portion of the 
spectra which was of interest. 

Continuous time series records of al] measurements were taken for 
roughly 2-1/2 hours each day. Record lengths of 24 minutes from each 
data set were subjected to analysis. Criteria which determined record 
lengths included economy of computer usage and resolution over the fre- 
quency range of interest. 

A maximum lag time was chosen as 5% of the record length. This re- 
sulted in a spectral band width resolution of 0.007 Hz and each spectral 
estimate having 40 degrees of freedom. The ninety percent confidence 
limits for 40 degrees of freedom using a chi-square distribution are found 
to be between 0.72 and 1.5] of the measured power spectral estimates. 

A mean value was computed for all data sets and the data was linearly 
detrended to remove tidal effects. Variance, standard deviation and 
average period were calculated. The average period was determined by 
calculating the time between zero upcrossings. The calculated average 
periods are lower than the visually observed periods due to a number of 
perturbations, such as noise and capillary waves, increasing the zero 
upcrossing occurrences. A probability density function for each data set 
was calculated and plotted. Comparisons were made with Gaussian and Gram- 


Charlier distributions using the chi-square goodness-of-fit test. Variance, 


2] 





standard deviation, skewness and kurtosis of the distributions were 
computed and are summarized in Table II. 

For each data set an auto-covariance function was determined and then 
smoothed with a Parzen window. Applying a Fourier transform to the smoothed 
auto-covariance function resulted in the power spectrum. An examination of 
the power spectrum shows the regions of greatest potential or kinetic 
energies and their respective frequencies. 

Cross-spectra were computed by calculating a cross-covariance between 
data sets, smoothing it with a Parzen window and applying a Fourier trans- 
form. Coherence and phase were then determined from the cross-spectrum. 

A coherence and phase versus frequency plot will indicate the regions 


and degree of linear relationship and phase between two data sets. 
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RESULTS 


A. QUALITATIVE DESCRIPTION 

Observation of various breaker types which occur on different beaches 
has resulted in a number of similarities being noted. A typical analog 
record of waves and onshore velocities obtained at Torrey Pines Beach is 
shown in Figure 4. In general, a sawtooth shaped profile is created when 
there is a quick drawdown of water just before the breaker arrives, fol- 
lowed by a steep, vertical leading edge, and a sloping profile toward 
the trailing edge. On the trailing edge, secondary waves are often ob- 
served. Secondary waves are harmonics of the primary wave frequency and 
denote strongly nonlinear waves. At Torrey Pines Beach the breaking 
waves are generally of the spilling variety which spill rapidly at the 


crest and move down the wave. 


B. MEAN VALUES 

A mean value was computed for all data sets and linearly detrended 
to remove tidal effects. The sea surface elevation mean values ranged 
from 0.67 to 1.33. Wave heights increased till the 18th and then de- 
creased during the remaining period data was taken, reflecting well the 
actual environment. 

Analysis of the mean values for both the onshore and alongshore 
components of flow was complicated by an apparent rip current. A rip 
current was frequently observed located just south of the array of wave 
gauges and flow meters. The on-offshore flow for the five days under 


consideration at Flow Meter #1 was always directed offshore at about 0.1] 
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m/sec. At Flow Meter #2 the horizontal water particle velocity was di- 
rected onshore at roughly 0.25 m/sec. The flow at Flow Meter #3 was 

away from the coast on 16, 17, and 18 March at slow speeds and towards 
the shore on 21 and 23 March moving very slowly. Concurrently, the long- 
Shore flow at Flow Meters #1, #2 and #3 was generally moving downcoast 
with the peak velocity varying from flow meter to flow meter from day to 
day. Similar results were reported by Huntley and Bowen (1974) for a 


nearshore circulation cell which they attributed to a rip current. 


C. PROBABILITY DENSITY FUNCTIONS 
I. Gaussian and Gram-Charlier Frequency Distributions 
The computed probability density function was compared to the 
Gaussian and Gram-Charlier distributions and tested for comparison using 
the chi-square goodness-of-fit test. The closer the chi-square fit 
parameter is to zero, the better the fit (Table II). 


The Gaussian probability density function is given by 


_١1_ о (хем)2/20 (1)‏ > (×٭ام؟ 
ое т‏ 


where g is the standard deviation and u is the mean. The Gaussian 
distribution is completely described by the mean and variance, and has 
zero skewness and a kurtosis equal to three. This approximates the 

sea surface in deep water conditions. However, nonlinearities found in 
the surf zone introduce skewness and kurtosis values that deviate from 


Gaussian and result in a distribution more closely approximated by the 


2 





Gram-Charlier distribution. This is to be expected in that the Gram- 
Charlier pdf is calculated using the assigned parameters from the data, 
and includes moments higher than the second. The frequency distribution 
of Wave Gauge #1 on 23 March is shown in Figure 5. Appendix D contains 
distributions calculated for other data sets. 

Distributions for all flow velocities are either unimodal or have 
less pronounced secondary peaks compared to the distributions of wave 
heights. The distribution of the onshore flow of Flow Meter #1 on 23 
March is shown in Figure 6. As noted in Table II, this results in nearly 
equal fit parameters with the Gram-Charlier distribution being generally 
smaller and thus a better fit. Fit parameters for velocity components 
are much closer to zero than are the values calculated for wave components 
since the flow meters do not experience the surface irregularities to 
the same degree as do the wave gauges. 

2. Skewness and Kurtosis 

The sea surface distribution was found to be positively skewed. 
This means that positive values are large, but less frequent than negative 
values which are more frequent, but smaller, indicating a greater amount 
of time below the mean water level. In other words, the crests are steeper 
and more peaked, whereas the troughs are elongated and flatter which is 
now waves in shallow water theoretically appear. 

It should follow that the wave-induced particle velocities would 
have distribution similar to the waves. However, of the 34 cases listed 
in Table II, 17 had positive skewness and 17 had negative skewness with no 


readily apparent pattern. Skewness of the onshore flow at Flow Meter #1 
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was always positive. Alongshore flow at the same location had negative 
skewness in three of four cases. Onshore water particle velocity at 
Flow Meter #2 had negative skewness for all six days. The skewness of 
alongshore flow at the same point was positive in five of six cases. 

The onshore component of flow at Flow Meter #3 had a positive skewness 
for 16, 17, and 18 March and a negative skewness for the following three 
days on which data was taken. Skewness of the alongshore flow at Flow 
Meter #3 was positive for four of the six days. The reason for this 
disparity is not known. Similar anomalous results have been found by 
Thornton and Galvin (1975). 

Kurtosis indicates the peakedness of each parameter. Just prior 
to breaking the waves achieve the greatest degrees of peakedness. Visual 
observations indicated that most of the waves were breaking at or near 
Wave Gauge #1 on 17, 18, and 19 March; this correlates with their high 
kurtosis values on these days. Values of kurtosis for all flow veloci- 
ties again had no apparent pattern with values ranging from 2.6 to 8.6. 


The reason for this is not known. 


D. SPECTRAL ANALYSIS 
1. Sea Surface Elevation and Water Particle Velocity 
The power, coherence, and phase spectra of Wave Gauge #1 and the 
onshore component of the flow of Flow Meter #1] on 23 March are shown in 
Figure 7. The spectrum of the wave surface elevation, which is char- 
acteristic of all the records, show a narrow banded peak at a frequency 


of 0.07 Hz, corresponding to a period of 14.3 sec. Also evident are 
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peaks at 0.14 and 0.2] Hz which appear to be harmonics of the primary 
peak at 0.07 Hz. These harmonics appear to be physical as was observed 
from the strip chart earlier. The harmonics probably have some energy 
contribution due to the Fourier computational technique. Appendix E 
contains spectra calculated for all other wave and onshore flow data sets. 

The coherence values were high, ranging to greater than 0.9 in 
the maximum energy portion of the wave. Coherence begins to fall off 
at about 0.5 Hz except on the 17th and 18th of March, the two days on 
which short-crested waves were noted. On these days the coherence begins 
to decrease at about 0.3 Hz. The decrease in coherence is probably due 
to the relative increase in turbulence, the general decrease in energy 
level, and wave energy spreading. The decrease in coherence due to wave 
energy spreading is discussed in Section E. 

The phase angle, which according to linear theory should be zero, 
had an average value of less than 20 degrees over the highly coherent 
band of prominent wave energy. The phase angle tends to increase slightly 
over the coherent band and becomes random for the noncoherent region of 
the spectrum. 

In order to compare waves and flow velocity, the wave profile 
spectrum was converted to a theoretical velocity spectrum for comparison 
with the measured flow velocity spectrum using the linear transfer function 


such that 
(Е) = |н(6)| 45, (ғ (2) 


The transfer function is given from linear wave theory 
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H(f) = 2 rf _coshk(htz) (3) 


sinh kh 


where 

k is the wave number, 

h is the mean water depth, and 

z is the depth of flow meter below the mean depth (m) 
Quantitatively, the results showed that linear theory overestimated 
wave-induced velocity spectral components by about 50%; in previous 


] (1976) where the waves were long crested and 


studies by Thornton et 
arrived at near normal incidence, the opposite was true, in that linear 
theory underestimated the magnitude of the wave spectrum. Tne difference 
in results is probably due to wave directionality, since when the wave 
approaches at any angle other than normal the flow meter measures only a 
component of the actual flow. 
2. Onshore and Alongshore Flow 

The cross-spectra between the onshore and alongshore components 
of flow were calculated with the intent of finding the relation between these 
two parameters. Power, coherence, and phase spectra of the onshore and 
alongshore flows at Flow Meter #3 on 23 March are shown in Figure 8. 
The spectra, which are representative of all the records, shows very 
little coherence between velocity components. The phase angle is ran- 
dom over the noncoherent band. 

The reason for the lack of coherence between the two orthogonal 
flow components is believed to be due to wave directionality. As was 


seen earlier the onshore component of flow appears highly wave-induced; 
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concomitantly, the alongshore flow component is at near normal incidence 
to the wave direction and thus has only a small wave-induced contribution. 
The variance of the alongshore component f flow is significantly less than 
the on-offshore component except at zero frequency, which reflects the 


mean longshore current variation. 


E. TURBULENT VELOCITY VERSUS WAVE-INDUCED VELOCITY 
The total velocity can be separated into components of a mean, plus 

wave-induced, plus turbulent velocity 

)4( .'نا + لا + نا - نا 
The wave-induced and turbulent velocity spectral components are assumed‏ 
to be statistically independent. For unidirectional waves aligned with‏ 
the measured velocity component, the co-spectra between waves and velocity‏ 
is then given by‏ 

Supt f) = Supt f. (5) 
Assuming statistical independence, the horizontal velocity spectrum 
is obtained from 

Su CF) - Sut) i Sy(f). (6) 
A further assumption is made that the waves and wave-induced velocities 


are given by a constant parameter linear process where the coherence 


is equal to unity 





Substituting (5), (6), and (7) into the definition of coherence between 


the total horizontal velocity and waves results in 
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S 
y 2 (ғ) - 1 + 2 (f) E SU(f) 
U (f) ШІ (8) 


Increasing lack of coherence is due to an increasing high ratio of turbu- 
lence. The coherence values indicate the percent of total velocity 
which is associated with the wave velocity. 
Wave-induced velocity can be determined using (8) 
2 


500?) = Уу (705,09). (9) 


The wave-induced velocity spectral component calculated in this way 
should be less than the actual value, the reason for this being that 
the linear coherence between waves and velocities will always be under- 
predicted due to nonlinearities which are always large in a breaking 
wave and also as a result of directional spreading of the incident 
wave energy (Battjes (1974) ). 

The turbulent velocity spectrum is derived by subtracting the cal- 
culated velocity spectrum from the measured velocity spectrum, 


5 (1) = S (f) - sun). (10) 


" 
The ratio of the turbulent velocity intensity to the wave-induced 
velocities intensity is shown in Figure 9. The data used for collapsing 
and plunging breakers, denoted by squares, was measured during earlier 
studies by Thornton et al (1976). The collapsing and plunging breakers 
were measured on beaches where the waves were long crested and arrived 


at near normal incidence. The stars represent data collected during this 


experiment. A value of the velocity intensity is given by taking the 
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square root of the variance which is computed by integrating the 


velocity spectrum, Su); across all frequency bands 


52 = а (па. (11) 
The measurements show a small spread of values for the ratio of turbulent 
to wave-induced velocity intensities. The average ratio calculated is 
approximately 0.82. This would indicate that, to a first approximation, 
the wave-induced kinetic energy is approximately 1.2 times the turbulent 
Kinetic energy. Thus the velocities in the surf zone on the average are 
highly wave-induced. 

If waves are not unidirectional, and generally they are not, it is 
necessary to account for directional spreading. The wave angle striking 
the wave staff and the type of wave, long or short crested, must be con- 
sidered. In general, short crested waves predominated at Torrey Pines 
Beach. On the 17th and 18th of March they were particularly noticeable. 
The coherence between surface elevation and velocity components of short 
crested waves is less than that of long crested waves due to directional 
spreading. Using the wave train model developed by Yefimov and Khristoforov 
(1971) the ratio of coherence between waves and onshore flow for the short 
crested waves to that for long crested waves is given by 


2 
۰۰۰۰۰۰۱۷٦00 ے0‎ 45, (8/35 17 (12) 


een (1 + Sy (F)/Sy (01/0059 117 


The model describes long crested waves using a train of two-dimensional, 
planar random waves traveling in direction 0, . This corresponds to a 


bivariate spectrum given by 
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o. (f,01) = s (f)8(0-6.) , (13) 


where 6(0-0 ) is the Dirac delta function. 
Short crested waves appear as a train of three-dimensional waves 
whose average direction of spectral components is assumed to be of the form 


2 
Y, (6) = | (2/т) 5,(Ғ) соз (0-6.)., 
0 for| 0 - Gol> 7/2, (14) 


Neglecting turbulence in order to examine the effect of directionality, 
the ratio of the coherence between waves and wave-induced velocities for 


short crested waves to that of long crested waves becomes 


т.е.) . 1/3ctn20,) ! 
E o 5 
791,91) 1 


As an example, when 6. equals an angle of 30°, the ratio equals 0.9. 


Conversely, the influence of turbulence on the ratio of coherence 
between waves and velocities for short crested waves to that of long 
crested waves can be examined by disregarding the directionality and 


letting 9=0. The ratio is given by 


2 -1 
602 (105 0073500 0 7 
У TASAS A T 


A summary of the ratios considering the turbulence and/or direction- 
ality factors at various angles is tabulated in Table III. The turbulence 
factor, SA Sy), is assumed for comparison purposes to be equal to 
0.75, the average value for the measured data in Figure 9 as denoted by 


the squares. 
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The coherence between waves and velocities was used to calculate the 
wave-induced velocity (equation 11). This method under-predicts the actual 
wave-induced velocity because of nonlinearities and wave energy spreading. 
Since the turbulent velocity is computed by subtracting the calculated 
wave-induced velocity spectrum from the measured velocity spectrum 
(equation 12) the smaller wave-induced velocity value makes the turbulent 
velocity value larger. The average ratio of the turbulent velocity 
intensity to the wave-induced velocity intensity was calculated to be 
0.825. In order to correct for the directional spreading it is necessary 
to multiply by the ratio У тб? 0M un (7,61). As an example, the 
ratio of the coherences for an angle of 15° is equal to 0.891. Multiplying 
this value by the average ratio results in a turbulence to wave-induced 


ratio of 0.75, which is equal to the earlier measured mean value. Hence 
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over the range of collapsing to spilling breakers a reasonable value for 


Ота 4 is approximately 0.75. 


F. SATURATION REGION IN THE SPECTRUM OF BREAKING WAVES 

Waves on Torrey Pines Beach characteristically broke as spilling 
breakers with the crest sliding down the face of the wave. When the 
fluid particles at the free surface move forward at a greater speed 
than the wave speed, c, breaking occurs. Strong harmonics observed 
in the wave and velocity spectra (Figure 7) are indicative of energy 
beina transferred from low to hiaher frequencies. The energy is even- 
tually dissipated by viscosity at the highest frequencies. When the 
transfer of energy is not rapid enough to balance the increase in energy 
density of the waves during shoaling, breaking occurs. The waves are 
said to be "saturated" with energy. Hence, a region of saturation 
would be expected through which energy is transferred from the low to 
higher frequencies. During shoaling and breaking of waves on a beach, 
the waves become saturated at low frequencies first imposing a bound 
on the peak energy density. 

Spectra of waves at breaking, or just inshore of breaking, are 
Shown in Figures 10 and 11 plotted on a log-log scale. The spectral 
estimates have been normalized by dividing by the total variance in an 
effort to bring the high end of the spectra to a single line. The 
spectral estimates have been block-averaged over frequency bands such 
that the log-energy values are linearly distributed on the log-frequency 
axis. The confidence intervals decrease with increasing frequency 


because of the increasing number of frequency bands averaged over. The 
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confidence interval is noted on the figures. The slopes of the log-log 
spectra at high frequencies for Wave Gauges #1 and Wave Gauges #3 сап 

be approximated by a -5 slope. This is indicative of a saturation or 
equilibrium region for wind generated waves in deep water. This was a 
surprise, as a -3 slope was expected for the waves measured in shallow 
water, Thornton (1976). А possible explanation is that most of the waves 
measured were waves inside the initial break point and would be classified 
as spilling breakers or reformed waves. It may be that these types of 
waves do not reach saturation condition as found for the more intense 
plunging and collapsing breakers. 

The normalized velocity spectra are given in Figures 12 and 13. The 
slopes of the velocity spectra at higher frequencies at Flow Meter #1 and 
Flow Meter #3 are both closely approximated by a -3 slope. A -3 slope 
is expected for the equilibrium region in both deep and shallow water 


and is consistent with the wave slopes. 
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Figure 12. Onshore water particle velocity spectra at 
Flow Meter #1. 
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V. CONCLUSIONS 


The values of the pdf's of sea surface elevation were positively 
Skewed verifying the asymmetrical shape of the waves in the breaker Zone. 
On the other hand, the skewness of the velocity components were both 
positive and negative and presented no readily apparent pattern. 

The mean values of the onshore and alongshore components of flow indi- 
cated the presence of a nearshore circulation cell probably associated with 
the rip current which was frequently observed to the south of the instru- 
ments. 

The values of the velocity energy-density spectral components calculated 
from wave spectra using linear theory indicate a qualitative, but not 
quantitative, relationship. In general, linear theory overestimated the 
magnitude of the velocity spectra because of wave directionality. 

The high coherence between waves and onshore flow indicated that most 
of the motion in the breaking wave is wave-induced. The coherence between 
waves and onshore flow was used to separate the turbulence to wave-induced 
velocity components. The measurements showed that over the range of col- 
lapsing to spilling breakers a reasonable value for the ratio of turbulent 
to wave-induced velocity is approximately 0.75. 

The waves lead the onshore flow in spectral phase on the average 
by less than 20 degrees, implying that the "curling" crest of the wave 
arrives prior to maximum water particle velocity. 

The wave and velocity energy-density spectra at higher frequencies 


exhibited slopes of -5 and -3, respectively. This is indicative of deep 
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water wave conditions. A possible explanation for the surprising spectral 
Slope is that most of the waves measured were either spilling breakers or 
reformed waves. It may be that the measured waves did not attain satura- 


tion conditions during breaking, as do plunging and collapsing breakers. 
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APPENDIX C 
CALIBRATION FACTORS 


DATE INSTRUMENT CALA* CALM** 
MARCH 
5:21,23 WG #1 -1.154 .00474 
7,19 WG #1 -1.284 .00503 
18 WG #1 -1.453 · 00520 
ALL FM £1 
DATES 
ONSHORE . 009 . 00764 
ALONGSHORE - .006 . 00747 
FM #2 
ONSHORE SA - . 00911 
ALONGSHORE - .084 -.00818 
WG #3 299 . 00156 
FM #3 
ONSHORE 2007 -.00820 
ALONGSHORE .006 -00763 


*Calibration Additive Factor for wave heights is given in meters; 
CALA for flow is given in meters per second. 


**Calibration Multiplication Factor for wave heights is given in 
meters per bit; CALM for flows is given in meters per second per bit. 
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APPENDIX D 
PROBABILITY DENSITY FUNCTIONS VS. 
GRAM-CHARLIER DISTRIBUTIONS 
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